An inverse method of designing cylindrical cloaks 
without knowing coordinate transformation 
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An inverse way to define the parameters of ideal cylindrical cloaks is de- 
veloped, in which the interconnection between the parameters is revealed 
for the first time without knowing a specific coordinate transformation. The 
required parameters are derived in terms of the integral form of cloaking 
generators, which is very general and allows us to examine the significance of 
the parametric profiles. The validity of such inverse way and the invisibility 
characteristics are presented in full-wave numerical simulation of plane wave 
scattering by cloaked cylinders. © 2009 Optical Society of America 

OCIS codes: 

Designing invisibility cloaks has received increasing interest from both engineering and sci- 
entific communities. Coordinate transformation is commonly used to control electromagnetic 
fields and render an object invisible to electromagnetic radiation [1,2]. This approach is gen- 
eralized from the cloaking in terms of the conductivity [3] , and further applied to acoustics 
and electromagnetics [4-6], which provides the possibility of concealing not only passive but 
also active objects within the interior [7,8]. 

Given a virtual space with unit parameters, coordinate transformation performs the tran- 
sition from the non-scattering virtual space to the realistic physical space, which in turn does 
not interact with the incidence. It is due to the fact that coordinate transformation method 
relies on the invariance of Maxwell's equations throughout the spatial transformation. This 
method can thus be applied to other classical waves (e.g., pressure wave [9] and surface 
liquid wave [10]) excluding elastic wave [11] provided that the symmetry invariance of the 
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corresponding wave equation is maintained under the coordinate transformation. Therefore, 
traditional cloaks need to be anisotropic and inhomogeneous. The anisotropy and the in- 
homogeity of the cloak, as the deformation of the space, bend the wavefront around the 
cloaked object and enable waves to emerge on the other side along the propagation direction 
without any disturbance. Based on this transformation technique, cylindrical cloaks with 
circular [12], elliptical [13] and arbitrary cross-sections [14] have been studied, and practi- 
cal attemps to realize the cylindrical cloak have been made with promising experimental 
results in microwave [15] and optical [16] regimes. Three-dimensional spherical electromag- 
netic cloak has been studied analytically by Mie theory [7] for first-order transformation and 
Von Neumann's method [17] for higher-order transformation, respectively. The realization of 
electromagnetic spherical cloaks by multilayered isotropic coatings has been proposed [19]. 

Traditionally, in order to derive the anisotropic parameters for the previous cloaks, one 
has to know the coordinate transformation beforehand. On the contrary, in this letter, we 
present a novel methodology to determine the required cloaking parameters for cylindrical 
cloaks without knowing specific coordinate transformations. Once any parameter is given, 
the rest parameters can be derived in integral form associated with the introduced cloaking 
generator. In this method, it is not even necessary to know the exact profile of the first 
given parameter. The cloaking generator together with the boundary condition, in fact, 
replace the corresponding spatial transformation. Full-wave simulation results are provided 
for verification. This work is an important step forward in the search of desirable cylindrical 
cloaks via parametric profiles. 

For impedance matching purpose, both the relative permittivity and permeability tensors 
of the cylindrical cloak are assumed to be equal, i.e., e(r) = JZ{r) = £(r) = ( r (r)rr + 
C v {r)(p$ + C e (r)zz. 

First of all, we only know that a cylindrical cloak is designed by compressing the virtual 
region r' < b into the physical region a < r < b via an unknown prescribed transformation, 
where the prime corresponds to the virtual space. The variables of a and b are the inner and 
outer radius of the cylindrical cloak in the physical space, respectively. It is assumed that 
the spatial compression is only with respect to the radial direction, and thus the Jacobian 
matrix is diagonal though we still have no information of the specific form of the coordinate 
transformation. As a result, one can obtain 
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denote three principal stretches of the Jacobian matrix. Then three equations can be derived: 



Cr(r)C„(r) = l, Cr(r)Ur) = -, — = ^(r)C(O- (3) 
By manipulating Eq. (3) and eliminating the term ( v (r), we derive the differential equation 



regarding radial and transverse parameters 



d 



ry/Cr(r)( z (r) 



dr 

Integrating Eq. (4), one has 



rVWCJr)— 1 ~ = <z(r). (4) 
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where C is the integration constant. Due to the spatial compression (r' = when r = a), it 
can be seen in Eq. (3) that Cr( r )Cz( r ) = 0, resulting in C = 0. Another requirement (r' = b 
when r = b) leads to the normalization condition for the ^-component parameter 

b 

2r 1 ( z (r 1 )dr 1 , (6) 

which plays an important role in finding cloaking parameters. Here, we introduce cloaking 
generator g{r) proportional to ( z (r), i.e., g(r) = C ( z (r) where C is an arbitrary constant. 
Then ( z can be presented in the form 

2 J a riginjdn 

Radial and azimuthal parameters can be expressed as 



r 2 g(r) ' 2 firig(n)dn' 

Only after all parameters are determined can corresponding coordinate transformation for 
such cylindrical cloaks be found, which is in contrast to the existing design approaches of 
transformation based cloaks: 

Ja r i9(ri)dn 

One may consider some interesting situations. When two permittivities out of three are 
equal (i.e., ( r = or ( z = ( v , or £ r = ( z ), it can be shown that in each situation only the 
trivial cloak (r' = r) is possible. Thus, ideal cylindrical cloak has to be realized for three 
different permittivities: e r ^ ^ e z [15,18]. This is different from the case of spherical 
cloaks with only two different parameters [1,19]. 
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Table 1. Cloaking parameters under different power cloak generators for ( z (r) profile. The 
symbol n is an arbitrarily positive constant. We introduce T n (r) = (r — a)(a + r + nr)/(2 + 
3n + n 2 ) to simplify the expressions corresponding to the second power profile. 
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In <73 (r), it is simply Pendry's classic cylindrical cloak when n = 1. 




Fig. 1. (Color online) The total electric field Re[^] on the x-y plane when n = 1 and n = 20 
for the power cloak generators. The plane wave is propagating along x-axis and its electric 
field is polarized along z-axis. The working frequency is 3GHz. The inner rod is a PEC at 
the radius of a = A, and the outer radius is b = 2A. 

To demonstrate the implementation of the cloaking generator and the proposed inverse 
design method, we just consider and focus on three different profiles of the power generators 
as given in Table 1. Certainly, one may design other profiles for the cloaking generator (e.g., 
parabolic, periodic, exponential, and even sinh profiles), which is out of the scope of the 
current letter. 

When n = 0, g\ and g 2 power cloak generators are identical, and their ( z components are 
constant. As an numerical example, we present two sets of near-field patterns for the power 
generators in Table 1, i.e., n — 1 and n = 20. In Fig. 1, it is obvious that at n — 1, g± power 
cloak exhibits larger disturbance in the near-field pattern not only in the wavefront but also 
in the variation of magnitude, and both g 2 and g% outperform g\ under the same condition 
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(the meshing of n = 1 is the same for g±, g 2 , and g^) simulated in COMSOL Multiphysics. If 
the meshing can be extremely fine, there will be no disturbance for all, while it would require 
too much computer memory in simulation. This non-ideality of the simulation conditions can 
be considered as a non-ideal cloak consisting of a number of discrete cylindrical layers. Also, 
at extremely high order (e.g., n = 20), g\ outperforms g 2 and g^ since it results in smaller 
forward scattering. The field distribution oscillates severely near the outer radius in g 2 and 
#3 power cloaks (n = 20), implying that the meshing has to be finer in this region. It is noted 
that the field inside the cloak is more homogeneous at n = 20 because at larger ns the rays 
inside the cloaking shell will propagate very close to the outer radius. From the transform 
point of view, if the order of the power cloak generator becomes larger, more portions of 
the original virtual space is transformed into the region close to the outer boundary of the 
physical space. In theory, the cloaking performance remains while the field distribution inside 
the cloak is more squeezed into the outer radius. Thus, it needs much finer meshing close 
to the outer radius, and the cases of n = 20 all adopt the same fine mesh. However, we 
cannot make the mesh extremely fine due to the limitation of computer memory. Moreover, 
considering that the order n can be any arbitrary value, there may exist an optimal range of 
n in each power cloak in which the invisibility performance will be further improved when the 
optimization is applied to the disretized model [20], i.e., the calculated total scattering width 
becomes smaller than that of Pendry's classic one under the same computational criteria. In 
other cloak designs other than power cloaks, it is possible to find such optimal cases as well. 
Although these characteristics are of technical importance, in our letter, we just focus on 
the proposed inverse method which is useful to design cloak generators and determine the 
optimal cloak in practice. 

In this letter, we reported an inverse method to derive the cloaking parameters for cylin- 
drical cloaks which does not need to know the required coordinate transformation first. It 
provides us larger degree of freedom in the design of cylindrical cloaks in practice. Instead 
of considering the coordinate transformation, we can directly envisage the importance of the 
profile of just one parameter, based on which the rest can be analytically and uniquely deter- 
mined. The numerical result confirms the validity of the proposed approach, and also brings 
the community's attention to an alternative design method for various cylindrical cloaks. 

We thank Prof. Steven Johnson and Prof. John Joannopoulos for their stimulating com- 
ments and discussion throughout the preparation of this paper. 
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